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We measured the effects of periodic surface holes, created using a focused ion beam, on the 
phase diagram of the vortex matter in high-T c Bi2Sr2CaCu2 0s+(5 crystals. Differential magneto- 
optical measurements show that the irreversibility line is shifted to higher fields and temperatures, 
with respect to the pristine melting line. The irreversibility line displays weak field dependence 
between integer matching fields indicating multiple- flux-quanta pinning at holes. We find reduced 
equilibrium compressibility of the vortex matter at integer matching fields, which is strong evidence 
for the existence of thermodynamic Mott insulator phases. Shaking with a transverse ac field 
surprisingly reveals first-order melting that is not shifted with respect to the pristine melting line 
and that seems to occur within the Mott insulator regions. This melting is understood to be the 
first-order transition in the bulk of the crystal beneath the surface holes. The transition is visible 
at the surface, despite the reduced vortex compressibility in the top layer. 



PACS numbers: 74.72.Hs, 74. 25. Ha, 74.25. Bt 

I. INTRODUCTION 

Experimental and theoretical studies of high-T c ma- 
terials with correlated pinning centers have led to the 
discovery of many novel phases of vortex matter, nonex- 
istent in the pristine materials. These new phases arise 
from the complex interplay between intrinsic point dis- 
order, correlated disorder, vortex- vortex interaction, and 
temperature. Correlated disorder in high-T c crystals is 
often introduced by heavy ion irradiation along the crys- 
tallographic c-axis, which leads to a random distribution 
in the a-b planed The effects of such random cor- 
related disorder on the high-T c phase dia gram are rel- 
atively well- u nderstood , both theoreticall y 3 I 4 I 5 I £ 1Z1£I and 
experimental!} ^ * 11 * 12 *. Measurements of crystals with 
periodic correlated disorder are limited, since it is not 
yet known how to physically realize periodic correlated 
disorder in thick samples. Experimental studies of peri- 
odic disorder must therefore choose between two options: 
study of thin samples, or study of artificial pins located 
at the sample surface. Efforts have focused mainly on the 
study of thin superconducting film s with periodic pinning 
center J 13 ' 14 ' 15 ' 1 ^ 17 ' 18 ' 19 ' 20 ' 21 ' 22 '^. However, thin films do 
not necessarily retain the thermodynamic properties of 
the bulk crystal, due to enhanced point disorder. Com- 
parison of the resulting vortex phases to the thermody- 
namic phases of the pristine material is thus usually not 
possible. 

Theoretically, the thermodynamics of vortices in the 
presence of random correlated disorder have been stud- 
ied extensively. Nelson and Vinokui 3 mapped the pinned 
vortex matter onto a system of quantum 2D bosons. 
They predicted the Bose glass transition from the low- 
temperature Bose glass phase, in which vortices are local- 
ized, to a higher-temperature delocalized vortex phase. 



They also discussed the possibility of an incompressible 
Mott insulator (MI) phase, when the magnetic induc- 
tion of the sample B exactly equals the matching field 
= p(j)o, where p is the density of the columnar de- 
fects. Radzihovsk}^ extended this model to include ad- 
ditional phases for B > B§\ a low-temperature weak 
Bose glass phase, in which both vortices residing at pins 
and those at interstitial sites are localized, an interstitial 
liquid phase at intermediate temperature, in which inter- 
stitial vortices are free to move but those at pinning sites 
are still pinned, and a homogenous liquid phase at higher 
temperature, in which all vortices are delocalized. For 
periodic pinning centers, the various Bose glass phases 
may be modified 2 -^. For periodic surface holes, even such 
a modified description is expected to be valid only within 
some finite depth from the surface of the superconductor. 

Simulations of two-dimensional systems containing pe- 
riodic pinning centers that allow only single-vortex oc- 
cupancy demonstrate commensurate states at integer 
matching fields nB^, with permitted values of n depend- 
ing on the geometry of the pinning centers^!! Solutions 
of Ginzburg-Landau theory reveal additional commensu- 
rate states with multi-quanta vortices^ 2 ^ for more gen- 
eral sample and pinning center parameters. Different 
melting scenarios have been demonstrated for triangu- 
lar and kagome arrays at low matching fields 2 ^, and for 
square pinning arrays both at and in-between matching 
fields 29 . The square pinning array at the first matching 
field displays three phases: a low-temperature pinned 
solid with square geometry, an unpinned ("floating") 
solid with triangular geometry, and a high-temperature 
liquid that lacks long-range order. At higher commensu- 
rate matching fields, the floating solid phase is not found, 
but an intermediate phase with mobile interstitial vor- 
tices similar to the liquid is observed. Incommensurate 




fields display a pinned phase at low temperatures with 
extra vortices located at interstitial positions, a phase 
at intermediate temperatures in which some vortex mo- 
tion is present with both interstitials and pinned vor- 
tices participating, and a phase at higher temperatures, 
in which all vortices are mobile. The temperature at 
which mobility is observed for incommensurate fields is 
lower than the melting temperature at the commensurate 
matching fields^. For the triangular and kagome geome- 
tries, melting at the first matching field involves a low- 
temperature pinned solid and a high-temperature liquid 
only. Intermediate-temperature phases, in which some 
or all of the interstitial vortices are mobile, are observed 
at higher matching fields^. These melting transitions 
are expected to be most relevant to high-T c supercon- 
ductors^. 

Direct imaging experiments of low-T c thin films with 
artificial periodic disorder have shown that highly or- 
dered vortex states exist at integer nB^ and frac- 
tional (p/q)B ( f ) matching fields, with n, q and p inte- 
g er j30 | 3i [ -Q ue £ Q th e g e ordered vortex states, such films 
hav e demonst rated commensurate effects in critica l cur- 
rent) 13 * 14 * 15 * 16 * , magnetization^, magnet o - resist anc d 1 8 * 1 9 \ 
and magnetic susceptibility measurements^- 21 * 22 * 23 * . Pos- 
sible phases and phase transitions of the vortex matter 
have been inferred from these measurements. Enhanced 
flux creep rate for B > was thought to be evidence of 
a transition from an incompressible MI state to an inter- 
stitial liquid stated Shapiro steps in transport measure- 
ments were understood to be a result of the coexistence 
of vortices pinned to artificial pinning sites and mobile in- 
terstitial vortices 3 - 2 ^. The behavior of the critical current 
was interpreted as evidence of two depinning energies, 
corresponding to the upper boundaries of the weak Bose 
glass and interstitial liquid phases^. Onsets of nonzero 
real and imaginary parts of the magnetic susceptibility 
were tentatively identified as the lower and upper phase 
boundaries of an interstitial liquid phased In addition 
to these states, the possibility of a saturation number 
n s > 1, corresponding to n s vortices at each pinning site, 
leads to multiple- quant a pinned vorte x states, which have 
been observed in many samples^ 2 ^. 

There are fewer experimental data regarding the ther- 
modynamic phases of high-T c superconductors with pe- 
riodic artificial pinning centers. The critical current 
in YBa2Cu307 (YBCO) thin films exhibited integer 
commensurate effects^ over a large temperature range. 
Scanning Hall probe measurements indicated that trap- 
ping of ~ 15 flux quanta is possible for 2.5 /im- 
diameter holes close to T c in YBCO 34 . Thin crys- 
talline Bi 2 Sr 2 CaCu 2 8+ 5 (BSCCO) samples with fully- 
penetrating periodic holes exhibited integer 35 and ratio- 
nal 3 -^ matching effects, in magneto-resistance and trans- 
port measurements, respectively. Similar samples with 
surface holes also displayed matching effects in magneto- 
resistance and critical current! 3 -^. A single study using 
thick BSCCO samples with surface holes displayed inte- 
ger matching in local magnetization 38 . In these studies of 





FIG. 1: (a) SEM image of the cross section of the surface 
holes. Hole depth is ~ 1.4 /im. (b) Part of one of the two 
arrays patterned on the sample, imaged by SEM. The distance 
between holes is 0.9 /mi. (c) A DMO image of the sample at 
T — 80 K and H — 21 Oe. The sample's edge is the vertical 
border indicated by the arrow. The two — 29.5 G arrays, 
Ml and M2 of 170 x 170 /xm 2 , appear as darker regions, due 
to their enhanced irreversibility. The irregular shape of Ml 
is a result of surface damage that occurred after writing the 
holes. 

BSCCO, the matching effects were visible in the field and 
temperature ranges at which the vortex matter is known 
to be in a liquid state in pristine crystals. No first-order 
melting step 39 was measured in these samples, and a full 
description of thermodynamic phases and transitions is 
lacking. 

In this work, we present an investigation of a thick 
BSCCO crystal, partially patterned with periodic sur- 
face holes created by a focused ion be am, mea sured us- 
ing differential magneto-optics (DMO) 40 41 42 accompa- 
nied by shaking with transverse ac fields We observe 
step-like behavior of the irreversibility line (IL), which 
may be a result of multi-quanta pinning to holes. We see 
a reduction in the DMO signal at integer matching fields, 
evidence of MI phases. We find a first-order melting tran- 
sition (FOT) in the patterned regions that is not shifted 
with respect to the pristine melting line. This FOT is ob- 
served even at integer matching fields, where the vortex 
matter in the surface layer is essentially incompressible. 
We believe this FOT to be the melting transition of the 
bulk of the crystal beneath the patterned surface. 



II. EXPERIMENTAL DETAILS 

Several samples were prepared and studied. Here we 
present a detailed investigation of a 2750 x 740 x 30 /im 3 
BSCCO crystal (T c ~ 90.5 K), with two triangular arrays 
of periodic holes patterned on the top surface using an 
FEI Strata 400 focused ion beam system. Figures [TJa) 
andQb) show SEM images of the hole profile and period- 
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icity, respectively. The measured hole depth was approx- 
imately 1.4 /im. Hole diameter decreases from ~ 0.6 /im 
at the sample surface to ~ 0.3 /im at a depth of 0.7 jam. 
The lattice constant of both arrays was 0.9 /im, corre- 
sponding to a matching field of = 29.5 G. The dimen- 
sions of each array were approximately 170 x 170 /im 2 . 

DMO measurements were performed by modulating 
the applied field iJ||c-axis by AH = 1 Oe while sweeping 
temperature T at constant H or scanning H at constant 
T. Each measurement point required averaging over k 
CCD camera exposures, first at H + AH /2 and then at 
H — AH/ 2, and calculating a difference image. Each 
DMO image is the average of m such difference images. 
Using fc, m ~ 10 with a typical exposure time of 0.3 sec 
yielded a typical modulation frequency of ~ 0.33 Hz. 
Values of dB/dH were derived from the DMO images 
by dividing the local light intensity by the intensity of 
some region far from the sample, where it was assumed 
that dB/dH = 1 G/Oe. For quantitative data analy- 
sis, intensities were spatially averaged over typical area 
of ~ 50 x 50 /im 2 . As described previousl y 40 ' 41 1 42 1 , the 
DMO measurement with field modulation is essentially 
equivalent to the measurement of the real component of 
the low-frequency local ac susceptibility as obtained e.g. 
by Hall sensors^. Figure [ljc) shows a DMO image of 
part of the sample, taken at T = 80 K and H = 21 Oe. 
The average brightness of the patterned areas Ml and 
M2 is lower than that of the neighboring pristine sample. 
This is due to an elevated irreversibility line (IL) in the 
patterned areas as described below. 



III. RESULTS 

We first inspect the IL of the patterned regions. The 
IL is important in the context of possible Bose glass 
phases because it is thought to be the dynamic mani- 
festation of the thermodynamic Bose glass transitiorP. 
In DMO measurements, reversibility of the vortex mat- 
ter is quantified by modulating the applied field by AH, 
and measuring dB/dH, the change in the local magnetic 
induction due to the modulation. Strong pinning results 
in dB/dH = 0, whereas full reversibility corresponds to 
dB/dH = 1 G/Oe. The irreversibility threshold in the 
following was chosen arbitrarily at dB/dH = 0.8 G/Oe, 
with Til (Hil) denoting the temperature (external field) 
at which this threshold is reached. We emphasize that 
the resulting IL reflects the response of the vortex system 
at low frequencies. Transport measurements or DMO 
with current modulation could possibly map out addi- 
tional boundary lines similar to the derealization line of 
vortices from columnar defects^, however, such measure- 
ments are beyond the scope of the present study. Figure [2] 
shows dB/dH for the pristine region (Fig. [2^ a)) and pat- 
terned region Ml (Figs. |2^b) and[2|c)), measured by T 
scans at constant H. 

Focusing on the pristine region (Fig. [2|a)), we see 
a series of sharp peaks in dB/dH with paramagnetic 
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FIG. 2: (Color online) dB/dH, the change in magnetic in- 
duction due to field modulation of AH = 1 Oe, for different 
sample regions, taken during T-scans. (a) A pristine region. 
The pristine melting line T m appears as a series of light gray 
(light pink) paramagnetic peaks in dB/dH with values above 
1 G/Oe. This particular run was carried out on a sparse grid 
in T and H of 0.4 K and 0.5 Oe (denoted by the lines in 
the lower left corner). The apparent discontinuity in T m is 
an artifact resulting from grid spacing being larger than the 
width of the melting peak, (b) dB/dH of region Ml in the 
vicinity of the first matching field, — 29.5 G. dB/dH val- 
ues are reduced compared to (a) and a narrow dip appears at 
H = 33 Oe. The location of the pristine melting line T m is 
denoted by black points, (c) dB/dH in the patterned region 
Ml over a wide range of T and H. Matching effects (denoted 
by arrows) are visible at H = 33, 64, 95, 126, and 157 Oe, con- 
sistent with integer multiples of the predicted = 29.5 G. 
The pristine melting line T m is denoted by black points. The 
patches in the data are a result of slightly differing setup pa- 
rameters for the different experimental runs. 



dB/dH > 1 (light gray, light pink online), correspond- 
ing to the first-order melting transition T m from a low- 
temperature vortex solid to a high-temperature vortex 
liquicPSEU. The black dots in Fig.j2lb) show the pristine 
melting line T m extracted from Fig.E[a). The patterned 
region Ml (Fig. [2jb)), in comparison, shows no FOT. It 
does, however, exhibit two notable features. First, the 
IL of region Ml is shifted to higher temperature. This 
is seen by focusing on Til (white contour) in Figs. [2^a) 
and|2|b). Til of the patterned region Ml is significantly 
greater than Tjl of the pristine region, which occurs 
at temperatures lower than the pristine T m . The sec- 
ond notable feature in Fig [2^b) is a narrow finger near 
H = 33 Oe, approximately 1 Oe wide, for 82 < T < 87 K, 
in which dB/dH of Ml is suppressed. 

Figure [2|c) shows dB/dH of patterned region Ml over 
a larger range of H and T. The pristine melting line 
T m is plotted as black dots for comparison. The IL of 
Ml is clearly shifted to higher temperatures. Sharp fin- 
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FIG. 3: (Color online) dB/dH for B < B^, measured during 
T-scans. (a) dB/dH of a pristine region. The melting transi- 
tion Tm appears as a line with dB / dH > 1 G/ Oe (light gray, 
light pink online). The temperature Til at which the IL is 
located (dB/dH = 0.8 G/Oe, white) is found below T m . (b) 
and (c) show dB/dH of patterned regions Ml and M2, respec- 
tively. Til of Ml and M2 is shifted to higher temperatures 
relative to the pristine Til- A sharp finger in Til of both 
Ml and M2 appears at H = 33 Oe, where B = B<f>. Negative 
values of dB / dH (black) correspond to negative permeability, 
due to geometrical barriers. 



gers, or narrow regions of H in which dB/dH of Ml 
is highly suppressed, occur at H = 33,64,95,126, and 
157 Oe (denoted by arrows), consistent with integer mul- 
tiples B/B^ = 1,2,3,4, and 5 of the predicted matching 
field B^ = 29.5 G. The minima in dB/dH as a function 
of H at matching fields are of both dynamic and ther- 
modynamic origin. Enhanced pinning at matching fields 
suppresses vortex motion, and hence also dB/dH. As dis- 
cussed below, these minima also indicate narrow ranges 
of H with reduced equilibrium compressibility, since the 
compressional modulus en is proportional to dH/dB^*. 
In addition to the fingers observed at integer B/B^, we 
see that Til between matching fields is step- like, with 
Til(H) weakly dependent on H between matching fields, 
and shifts in Til(H) occurring at matching fields. For ex- 
ample, Til ~ 76.5 K for 3 < H/B^ < 4, and T IL « 74 K 
for 4 < H/Bcf, < 5 (white contour). 

We now focus on the IL of the patterned regions be- 
low and in the vicinity of the first matching field. Fig- 
ure [3] shows dB/dH measured simultaneously in the pris- 
tine region and in the patterned regions Ml and M2, for 
B < B4. Til (dB/dH = 0.8 G/Oe, white contour) of 
the pristine region is located below the pristine melt- 
ing line T m . In region Ml (Fig. [3fb)), T IL is shifted 
to higher temperatures at all fields when compared to 
the Til of the pristine region (Fig. j3^a)). An additional 
sharp finger in the IL, extending to Til — 87 K, occurs 
at H = 33 Oe, or B = B^. The IL of patterned re- 
gion M2 (Fig. [3^c)) is shifted somewhat less, yet it too 
displays a sharp shift to higher T at the matching field 
B^. The reason for the difference in the shift of the IL 
of the two arrays may be a result of the difference in 
the arrays' locations (M2 is closer to the sample's edge) 
or due to some difference in the holes of the two arrays, 
which are not identical and may have different pinning 
properties. Still, for both arrays the IL is shifted upward, 



with an additional sharp finger at B^. In BSCCO crys- 
tals irradiated with low concentration of CDs the sharp 
finger is absent. Instead, a kink is observed in the vicin- 
ity of B^ which is believed to be the result of depinning 
of two different vortex populations. Below B^, Til is 
the temperature at which vortices located at CDs depin. 
Above £>0, Til is the depinning temperature of intersti- 
tial vortices^. Also, fractional matching features have 
been observed 37 in BSCCO samples with periodic sur- 
face defects. We do not detect any fractional matching 
features in the IL of the patterned regions; the reason 
for this is not clear, but is consistent with results shown 
elsewhere^. 

The IL in BSCCO is kno wn to be a dynamic feature 
of the phase diagram^ 4 * 47 * 48 *. However, since it indicates 
a region of the phase diagram in which there is a change 
in the system's dynamic response, it may indicate an un- 
derlying thermodynamic transition that occurs at simi- 
lar values of field and temperature. In order to detect 
a possible underlying FOT, and to determine whether 
the reduced dB/dH regions are truly a thermodynamic 
feature of the phase diagram, we studied the behavior 
of the first matching field in the presence of shaking. 
The shaking technique^ utilizes an in-plane, ac mag- 
netic field. It is known to suppress hysteretic behavior 
in BSCC O, enab ling the observation of thermodynamic 
prop ert ie J 48 | 49 | 5Q J . 

The effect of shaking is shown in Fig. [4j for the pris- 
tine and patterned regions Ml and M2 at 77 K. Increas- 
ing shaking amplitude H^ c leads to a systematic increase 
in dB/dH, which saturates near 1 G/Oe, as expected 
for a fully penetrable sample. A feature common to all 
three plots in Fig. [4j and thus unrelated to the surface 
holes, is an abrupt change from zero to negative dB/dH 
at ^ 8 Oe. Below 8 Oe, the sample is in the Meiss- 
ner phase, with B = 0. dB/dH < immediately above 
the Meissner phase corresponds to negative local perme- 
ability 51 . This effect, which occurs in BSCCO samples 
with platelet geometry, is a result of the geometrical bar- 
r i e j52 | 53 | anc j modulation of the vortex dome during 
the modulation cycle of the applied field H±AH/2. This 
negative permeability is also visible as a black strip at low 
fields in Fig. |3ja). It is interesting to note that the neg- 
ative dB/dH values are not visible at the highest shak- 
ing amplitude, H^ c = 121.8 Oe. This indicates that the 
shaking field enables the vortices to overcome the geo- 
metrical barrier. Consequently, the local negative per- 
meability changes to high positive local permeability, as 
seen in Figs.[4^a) and|4^c). This is the expected behavior 
in the absence of geometrical barriers, as demonstrated 
for prism-shaped samples^. 

There are two notable differences between the pristine 
(Fig. Qa)) and patterned (Figs. Qb) and[4^c)) regions. 
The first difference is the appearance of the matching fea- 
ture in the form of a dip in dB/dH near 32 Oe, denoted 
by arrows in Fig. [4j Shaking extends the range of tem- 
peratures for which this dip is visible well into the vortex 
solid region below T m . Without shaking, the first match- 
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FIG. 4: (Color online) dB/dH vs. applied field H (scanned up 
and down) for different values of in-plane shaking amplitude 
Hac, at T = 77 K for the pristine region (a) and patterned 
regions Ml (b) and M2 (c). Arrows denote the first matching 
field. Data are shown for different values of applied H^ c (from 
bottom to top): (>), 20.3 (V), 40.6 (A), 81.2 (0), and 121.8 
(□) Oe. Shaking frequency was 15 Hz for all measurements. 
Symbols appear every 37 data points. 



ing feature is not visible at this temperature (77 K, see 
Fig. [2|, due to the enhanced pinning in the vortex solid. 
With increased H^., the matching feature appears first as 
a step (Fig.^b), H^ c = 20.3 and 40.6 Oe) and then as a 
dip in dB/dH (H^ c = 81.2 Oe). In some cases, as shown 
in Fig. |ijc) for high H^ c , a peak appears in dB/dH im- 
mediately before the dip. The slight downward shift in 
H of the dip for increasing H^ c probably results from the 
increased penetration of magnetic induction B at higher 
H^c resulting in the same B = B^ at slightly lower val- 
ues of applied field H. The dip in dB/dH corresponds 
to a reduction in the compressibility of the vortex matter 
at B^. 

The second difference between the pristine and pat- 
terned regions in Fig. [4] can be seen away from the 
matching field. For the pristine region, the values of 
dB/dH increase gradually from dB/dH ^ at low field 
to dB/dH = 1 G/Oe at sufficiently high applied field 
H. For the patterned regions, the behavior of dB/dH is 
plateau-like, with the matching feature dividing between 
neighboring plateaus. This can be seen in Fig. gb) for 
H^ c = 81.2, 121.8 Oe, and in Fig.gc) for H^ c = 40.6 Oe. 
These plateaus are consistent with the observed step-like 
behavior of Til during T-scans, as shown in Fig. [5|c). 
The plateaus in dB/dH appear between matching fields, 
where Til is almost independent of H. The step between 
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FIG. 5: (Color online) dB/dH of region M2 as a function of 
H for different temperatures T and shaking amplitudes H^ c . 
The effect of shaking is similar to the effect of temperature 
(see text), (a) Matching effects at T = 74 K, B/B^ = 2,3, 
and (b) at T = 70 K, B/B^ = 3,4. H^ c = (V), 20.3 (A), 
40.6 «>), and 81.2 (□) Oe. (c) Matching effects at T = 74 
(A), 76 (0) and 78 (□) K, B/B^ = 2,3, without shaking, 
(d) Matching effects at T = 68 (>), 70 (V), 72 (A), 74 (0), 
and 76 (□) K, B/B^ = 3,4, without shaking. Arrows denote 
matching fields. Symbols appear every 30 data points. 



the plateaus appears at H = B ( 
steps in Til that occur at H = nB, 
for 



consistent with the 
The value of dB/dH 
each plateau in Figs. 4^b) andEFc) increases with in- 
creasing H^c- Figures [5 a) and[5[bj show the effects of 
shaking with different values of H^ c for B/B^ = 2, 3 and 
B/B^ = 3,4, respectively, for patterned region M2. The 
data for Ml are similar. Clearly the same matching fea- 
tures that appear for B/B^ = 1 , namely a step in dB / dH 
for low H^c that develops into a dip for higher H^ c1 are 
visible also for higher matching fields. 

Interestingly, increasing T and increasing shaking am- 
plitude H^c have similar effects on dB/dH immediately 
below the IL. This can be seen by comparing Figs. [5|a) 
and [5|c), in which dB/dH of region M2 in the vicin- 
ity of B /Bff, = 2, 3 is plotted for different values of H^ c 
and T, respectively. Increasing H^ c and increasing T 
(bottom to top curves) both tend to increase dB/dH, 
and both have the effect of transforming the matching 
feature from a step to a dip. However, the overshoot 
in dB/dH immediately below the matching feature ap- 
pears only at nonzero H^.. Similar behavior is observed 
for B/B^ = 3,4 in Figs.jtjb) and^d). 

We now address the question of first-order melting 
within the patterned regions in the presence of shak- 
ing. For the results shown below, we applied a 15 Hz, 
H^c = 81-2 Oe shaking field. We find that shaking shifts 
the IL to lower fields and temperatures and thus enables 
the observation of a FOT in the patterned regions of the 
sample. Figure [6] shows detailed scans of the H—T region 
in which the pristine T m line intersects the B/B^ = 1 
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FIG. 6: (Color online) dB/dH as a function of H and T near 
the intersection of melting and B/B^ — 1 matching, without 
(top panels) and with (bottom panels) shaking. Shaking pa- 
rameters were 15 Hz and 81.2 Oe. Results are shown for the 
pristine sample ((a),(b)) and patterned regions Ml and M2 
((c), (d) and (e),(f), respectively). 
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FIG. 7: (Color online) dB/dH as a function of H and T near 
the intersection of melting and B/B^ — 2 matching, without 
(top panels) and with (bottom panels) shaking. Shaking pa- 
rameters were 15 Hz and 81.2 Oe. Results are shown for the 
pristine sample ((a),(b)) and patterned regions Ml and M2 
((c), (d) and (e),(f), respectively). The slight curving of the 
matching effect to lower H for higher T in (c),(d),(e) and (f) 
is due to the increased penetration B(H) at higher T. 



IV. DISCUSSION 



matching line. dB/dH is shown for both patterned re- 
gions and for the pristine region, without and with shak- 
ing (top and bottom panels, respectively). For the pris- 
tine region (Figs.[6ja) and [6|b)), dB/dH at lower T and 
H is raised slightly by shaking, and the pristine melt- 
ing line T m , which appears as a line with paramagnetic 
dB/dH > 1 G/Oe (light gray, light pink online), remains 
essentially unchanged. For the patterned region Ml, no 
FOT was visible without shaking (Fig.[6jc)). With shak- 
ing (Fig. |6|d)), a FOT became visible. It appears to be 
located at the same temperatures and fields as the pris- 
tine melting line T m . Remarkably, the T m line is clearly 
visible even at the bottom of the B^ matching dip. For 
the patterned region M2, shaking was not needed to un- 
cover the FOT (Fig.^e)). While shaking raised dB/dH 
values overall, it did not change the location or the na- 
ture of the FOT (Fig.gf)). Similar results are shown for 
B/B^ = 2 in Fig. [7] In this case, shaking was necessary 
to view the FOT in both patterned regions Ml (Fig.[7jc) 
andfjjd)) and M2 (Fig.^e) andgf)). Note that the lo- 
cation of the FOT of the patterned regions in the phase 
diagram is indistinguishable from the location of the pris- 
tine melting line, T m . Moreover, T m and the nB^ lines 
seem to intersect with no apparent interaction, as if the 
periodic pinning potential of the holes has no effect on 
melting. No additional FOT was detected for either of 
the patterned regions. We emphasize that at the points 
in the phase diagram where the FOT meets the matching 
fields, a contradictory behavior of the vortex lattice oc- 
curs. On one hand, at the FOT there is a jump in vortex 
density. On the other hand, at matching fields the vortex 
matter exhibits a strongly enhanced compressibility mod- 
ulus en ~ (dB / dH)~ x that exists both below and above 
T m . This apparent contradiction is discussed below. 



In order to understand the observed behavior of the IL, 
we consider two possible physical scenarios 55 . In the first 
scenario, shown schematically in Fig. [8^ a), we assume 
that each hole can pin only a single vortex. As a result, 
two vortex populations are present for B > B^\ vortices 
located at holes, and interstitial vortices located between 
holes. The interstitial vortices are subject to a caging po- 
tential caused by the vortices located at holes^, that is 
assumed to be weaker than the pinning potential at holes, 
but stronger than the pristine pinning. This gives rise to 
a depinning transition of the interstitials, which we iden- 
tify with T/l, at a temperature above the pristine melting 
temperature T m . Alternatively, we consider a scenario in 
which there is multi-quanta pinning by holes. We assume 
that below the IL, all vortices are located at holes, while 
above the IL, some vortices are depinned from holes, and 
thus mobile, as shown schematically in Fig.[8|b). Due to 
repulsion between pinned vortices, the pinning force per 
vortex is expected to decrease as a function of the num- 
ber of vortices pinned to the hole 5 -^. We therefore assume 
that the pinned vortices residing at holes depin one at a 
time, as T is increased. Within this multi-quanta sce- 
nario, Tjl corresponds to the temperature at which the 
first vortices depin from holes. Either of the two scenar- 
ios must provide an explanation for the observed behav- 
ior of the IL: the plateaus in the IL away from matching, 
the shift in the IL to lower T and H in the presence of 
shaking, and the sharp dips at the matching fields. 

We begin with the plateaus and steps in the IL. In 
the multi-quanta scenario, the maximum number of vor- 
tices per hole n max (T) is determined by temperature- 
dependent hole pinning strength and repulsive interac- 
tions between vortices located at the hole, and decreases 
with increasing temperature 5 -^. We denote the tempera- 
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FIG. 8: (Color online) (a) and (b): A schematic view of the 
vortex matter for B/B^ = 3, within (a) the single- vortex pin- 
ning scenario, and (b) the multi-quanta scenario. Gray circles 
indicate holes. Black dots indicate vortices, (a) In each unit 
cell of the pinning lattice, a single vortex is pinned to the hole 
and two vortices are interstitial. Below (above) the IL intersti- 
tials are pinned (mobile), (b) In each unit cell, two vortices re- 
main pinned to the hole. Below the IL the third vortex is also 
pinned to the hole (not shown). Above the IL, it is depinned 
and mobile, (c) and (d): A schematic description of the steps 
in the IL, within the multi-quanta scenario, in the absence (c) 
and presence (d) of thermal fluctuations or shaking, (c) The 
maximum number of quanta per hole n max (T) (solid black 
line) is expected to decrease as a function of T, resulting in 
temperatures T n (dashed lines) above which holes pin n — 1 
vortices only, and additional vortices become mobile intersti- 
tials. Thus for each interval (n — 1) < B /B^ < n, Til =T n . 
(d) Thermal fluctuations shift Til (solid line) from T n (dashed 
lines) to T n < T n . At B/B^ = n, the vortex lattice remains 
stable to the temperature T n due to a thermodynamic MI 
phase. Resulting phases include a multi-quanta-solid (MSQ) 
below Til, an interstitial liquid (ILq, white) above the IL that 
terminates at Ti, a homogeneous liquid (HL) above Ti, and 
MI "fingers" at B/B^ = n (black). 



ture at which n max decreases from n to n — 1 by T n (see 
dotted lines in Fig. |8jc)). We assume that for the ap- 
plied fields H < 6£?0, all vortices are pinned to holes at 
sufficiently low T. As T is increased above T n , holes may 
only pin n — 1 vortices. Therefore vortices abruptly depin 
from holes occupied by n vortices, leaving n — 1 vortices 
per hole. The depinned vortices are mobile, leading to a 
fast onset of reversibility. We therefore identify Til = T n 
for (n — 1) < B/B^ < n. The resulting IL thus displays 
steps and plateaus, as shown schematically in Fig. 8|c). 



Finite temperature and slight hole variability are likely to 
cause some variation in n max (T) at different holes. This 
would lead to some smearing of the irreversibility tran- 
sition and to a weak dependence of Til on B between 
matching fields due to different mixing of the T n as field 
is varied. This schematic description neglects the effects 
of thermal fluctuations, which will be discussed later on. 
The IL shown in Fig. [2|c) (white contour) indicates that 
Til displays three discrete steps in the temperature range 
T = 72 — 77 K. The weak temperature dependence of 
the IL between matching fields, as well as the plateaus 
in dB/dH in Figs. [4] and [5j seem to indicate that the 
IL is not strongly affected by interactions between vor- 
tices at neighboring holes. Rather, it is governed by the 
hole pinning energy and the repulsion between vortices 
within a single hole, leading to the step-like Til- Within 
the single- vortex pinning scenario, in contrast, Til(H) 
is the depinning temperature of interstitial vortices, that 
is expected to decrease rather smoothly with field as the 
density of interstitials increases. 

We now address the the shift of the IL to lower temper- 
atures in the presence of shaking. Shaking and increased 
thermal fluctuations seem to have a similar effect on the 
IL (see Fig. [5J. Within the multi-quanta scenario both 
provide a mechanism for hopping of vortices between va- 
cancies at holes, thus increasing the dynamic dB/dH and 
decreasing Til- Between matching fields and slightly be- 
low Til, the holes are below their full pinning capacity, 
so shaking may provide a way for vortices to hop between 
holes more easily, or to depin from holes and move to in- 
terstitial positions. The plateau-like behavior indicates 
that this shaking-induced hopping is roughly indepen- 
dent of the number of vacancies that exists. Instead, 
the degree of hopping or depinning is dependent on the 
balance between the roughly field-independent pinning 
energy of the n max th vortex and the activation energy 
of the applied shaking. Within the single- vortex pinning 
scenario, above the IL the interstitials are mobile, there- 
fore immediately below the IL they are weakly pinned. 
Shaking will thus assist in overcoming the weak pinning 
potential, decreasing Til- The existence of the plateaus, 
however, cannot be easily explained in this case. 

We now address the third experimental finding, namely 
the finger-like dips in dB/dH at matching fields. Focus- 
ing on Figs. [6] and [7| we see that unlike the rest of the IL, 
the temperatures at which the fingers terminate are not 
affected by in-plane shaking. This strongly suggests that 
unlike the IL, which is a dynamic feature of the phase dia- 
gram, the fingers at matching fields are a thermodynamic 
feature. A thermodynamic minimum with dB/dH = 
would indicate a plateau in the equilibrium B(H) and 
a diverging bulk modulus en ~ dH/dB*^, clear signa- 
tures of the incompressible MI phase. The finite, positive 
minima observed in dB/dH (see Figs. [IJb), [4]^c), and[5| 
correspond to a reduction in the positive slope of B(H) 
at nBtf). These finite values at nB^ may be a consequence 
of the finite size of the patterned region, which prevents 
infinite divergence of en, or of the broadening of dB/dH 
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due to the modulation of AH = 1 Oe, which is at least 
as wide as the width of the dip. For samples with ran- 
dom defects, it has been argued that the MI phase is 
destroyed by repulsive vortex interactions, possibly re- 
taining "lock -in" effects such as a finite peak in the bulk 
modulus c n 46 57 . In our measurements, however, the 
pinning is ordered. At matching fields, pinning energy 
and vortex-vortex interactions both stabilize the vortex 
latt ice. I n this case, observation of a MI phase is pos- 
sibl e 146 1 571 . We believe that the sharp matching features 
we observed in dB/dH at nB^, that are not affected by 
shaking, are a strong indication of thermodynamic MI 
phases. 

The fingers of reduced dB/dH, or MI phases, may be 
understood in the context of the multi-quanta scenario. 
For (n — 1) < B/B^ < n and T < T n , in the absence 
of thermal fluctuations, all vortices are pinned to holes 
(see Fig. |8jc)). The maximum number of vortices that 
can be pinned nB^ / O is greater than the actual number 
of vortices B/(j)o, resulting in below- full occupancy, or 
"vacancies", at some of the multi-quanta holes. Thermal 
fluctuations (which were neglected in the schematic de- 
scription in Fig.^c)) are expected to lead to vortex hop- 
ping between the vacancies, resulting in enhanced vortex 
mobility. Thus, for (n — 1) < B/B^ < n, vortex dy- 
namics lead to a reduction in Til from T n (Fig. [8^d), 
dashed lines) to some lower temperature T n . Exactly 
at matching, an additional thermodynamic consideration 
enters. The total hole capacity equals the number of vor- 
tices, and therefore there is a finite energy cost for adding 
an extra interstitial vortex. As a result, the equilibrium 
B(H) will exhibit a plateau over a finite range of H, that 
is not affected by shaking. The corresponding minima in 
dB/dH, or equilibrium MI phases, may thus be observed 
up to T n (black fingers in Fig. |8jd)). 

The fingers of reduced dB/dH may also be partially 
understood within the single- vortex pinning scenario. At 
matching fields, we expect both the vortices at holes and 
the interstitials to be ordered in a configuration commen- 
surate with the hole array, that enhances the pinning 
potential and may result in a thermodynamic MI phase. 
Although this picture is correct for commensurate match- 
ing fields (B/B^ = 1,3,4,7, etc. for a triangular array 
of holePl), it is not strictly correct for incommensurate 
matching fields. At B/B^ = 2, for example, there are 
two interstitial positions with equal energy per unit cell. 
This would imply that hopping of interstitials is possi- 
ble also at incommensurate matching fields, resulting in 
greater interstitial mobility, and hence a weaker matching 
effect. This is not what we observe in Fig.|2|c), where the 
fingers appear similar at the first and second matching 
fields. 

We now consider the theoretical plausibility of the 
multi-quanta scenario. The condition given for multiple- 
quanta pinning in holes i^S r 3 > £A 2 for A <C d, where 
r is the h ole radius , d is the inter-hole distance, and 
A = Ao/\A — T/T c is the London penetration depth 4 . 
For r = 0.2 /im, and Aq = 0.15 //irP^, we may expect 



multiple-quanta vortices for T < 88 K for our sample. 

We now estimate the number of vortices pinned at a 
hole as a function of temperature. We consider two physi- 
cal possibilities, following RefPH The first is that the nth 
vortex is pinned as long as there is some force near the 
edge of the hole that pushes it inwards. This corresponds 
to a pinned state that may be metastable, depending on 
its free energy. The saturation number 

n s « r/(20 (1) 

is the maximum n for which this occurs. The second 
possibility for multiple-quanta pinning is to require an 
equilibrium pinned state, namely, the free energy of the 
nth vortex located in the hole must be lower than its free 
energy far from the hole. In this case, the number of 
pinned vortices is given b}J^ 

1 i r M 2A 

n ^-ln— ;/ln . (2) 

2 2£ 1 1.78r v J 

Both n s and no decrease with T, which is consistent with 
the observed downward steps in Hil(T). Substituting 
r = 0.2 /im, £o = 2 nm, Ao = 0.15 /mi, and T c = 90.5 K, 
we obtain n s = 1 — 21 for T = 90.46 K down to T = 
74.5 K and n = 1-2 for T = 85.2 K down to T = 76.1 K. 
Note that Eq. [2] was derived from the pinning energy 
of a single hole. In the case of an array of holes, one 
should compare the free energy of a pinned nth vortex to 
its free energy at the midpoint between two neighboring 
holes. At this midpoint, there are positive contributions 
to the free energy from the neighboring occupied holes. 
Thus a higher free energy of the pinned vortex, or n > 
no, would still be an equilibrium state of the system. 
Eq. [2] should therefore be considered a lower limit on the 
equilibrium number of vortices pinned to a hole. The 
fact that interactions should raise estimated occupation 
numbers was also noted by^. 

We now compare these theoretical estimates to the ex- 
perimental values of n™^. We observed three decreasing 
steps in T IL (see Fig.^c)), with f 6 ~ 72 K, f 5 ~ 75 K, 
and T4 ~ 77 K. According to the schematic phase dia- 
gram plotted in Fig.|i|d), the values of T n are lower than 
the temperatures T n at which the n-th multi-quanta vor- 
tex depins in the absence of thermal fluctuations. From 
the schematic description in Fig. [8^c), we see that the 
difference T n — T n may be estimated from the differ- 
ence in the Til at, and slightly away from, B/B^ = n. 
From Figs. [6jc) and[7Jc), we estimate this difference to 
be T n — T n ~ 2 K. Tq,T$, and T4 are thus 74,77 and 
79 K, respectively, or equivalently, n max (T < 74) = 6, 
(74 < T < 77) = 5, (77 < T < 79) = 4, 
and n max (T > 79) = 3. We find that the extracted 
nmax values are much lower than the estimated n s , and 
slightly higher than the estimated no- Thus the equilib- 
rium multi-quanta pinning scenario described by Eq. [2] 
is more plausible. Note, however, that Eqs. [l] and [2] 
are based on the assumption of fully penetrating holes; 
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the number of vortices trapped by surface holes may be 
lower^HU. 

Summarizing our discussion of the IL, the multi-quanta 
scenario provides a more consistent explanation as com- 
pared to single- vortex pinning, both for the steps in the 
IL and for the similar MI fingers at both commensurate 
and incommensurate matching fields. Indeed, recent sim- 
ulations of BSCCO with surface holes similar to those in 
the experiment^ indicate that multiple-quanta occupa- 
tion of surface holes does occur, and that depinning may 
occur directly from holes. 

Finally, we address the apparent contradiction of ob- 
serving a FOT within the patterned regions at matching 
fields B/B^ = 1,2, as shown in Figs. [6] and [7] The ob- 
served FOT occurs at the same field and temperature 
values as the pristine solid-liquid transition T m , and we 
therefore assume that both transitions are of similar na- 
ture. This, however, seems to contradict the existence of 
MI phases at matching fields, since in the MI phase the 
vortex lattice is ordered, pinned, and incompressible, up 
to temperatures well above the pristine T m . 

The observed FOT can be understood qualitatively by 
taking the bulk beneath the surface holes into account. 
The surface holes are only ~ 1.4 /am deep, whereas the 
sample is 30 /am thick. Although the depth at which 
vortices are still sensitive to surface patterning is not 
known exactly, magnetic decorations of BSCCO crystals 
with square Fe periodic surface patterns indicate that 
the hexagonal structure of the lattice is recovered fully 
just 4.5 fim beneath the pinning potential at low tem- 
peratures^. Thus one may expect that sufficiently deep 
below the upper surface, the vortex "tails" will undergo 
a first-order melting transition at the pristine T m . At 
matching fields, however, the "tips" of the vortices at 
the surface are pinned to the periodic surface holes, and 
therefore behave as an incompressible solid. The result- 
ing situation immediately above T m is rather unique: the 
vortex "tips" are in a solid MI state, while the vortex 
"tails" are liquid. At the FOT, the vortex density in the 
bulk increases by A£?/(/> , where AB is the typical step 
in B at T m . If the "tips" of the vortices were fully in- 
compressible, this AB in the bulk would be completely 
shielded and unobservable at the top surface. Our data 
indicate, however, that within our experimental range of 
parameters the compressibility is high, but finite. Hence 
the MI top layer is sufficiently transparent to allow ob- 
servation of the paramagnetic peak at the FOT of the 
underlying vortex "tails". We conclude that upon in- 
creasing temperature, the FOT observed within the pat- 
terned region at matching fields indicates a transition 
from a solid bulk with an incompressible solid surface to 
a unique state of an essentially incompressible solid crust 
of vortex "tips" concealing a vortex liquid in the bulk. 

It is interesting to note that this solid crust could be 
used as a tool to study the interesting possibility of sur- 
face melting that was suggested to exist in the vortex 
lattice®^, similar to surface wetting in atomic solids. 
This will require investigating crystals of different thick- 



nesses and samples patterned with holes at both the top 
and bottom surfaces. 

No additional FOTs were observed, even in the pres- 
ence of shaking. This indicates that the depinning line 
of the interstitial vortices, and the derealization line of 
the vortices pinned to holes, are apparently not FOTs. 
More accurate magnetization measurements in the pres- 
ence of shaking are needed to check for the existence of 
second-order thermodynamic transitions. 



V. CONCLUSIONS 

BSCCO crystals with arrays of surface holes were mea- 
sured using differential magneto-optics accompanied by 
an in-plane shaking field. We observe reduced dB/dH 
of the patterned regions at integer matching fields B = 
nBfi. These features are extremely narrow, with a width 
in H of less than 1 Oe. The region of reduced dB/dH 
extends up to 87 K for the first matching field, termi- 
nating 3.5 K below T c . Shaking allowed the equilibrium 
matching feature to be observed both above and below 
the pristine melting line of the sample. This observa- 
tion is in contrast to previous dynamic m easure ments of 
BSCCO crystals with periodic surface hole d 35 * 3 ^ , in which 
matching effects were observed only well above the pris- 
tine melting line. Our finding of a sharply suppressed 
compressibility of an equilibrated vortex lattice at the 
matching fields is a strong evidence of the existence of 
Mott insulator phases. 

We observe a first-order melting transition within the 
patterned areas, both away from matching and at the 
first and second matching fields. Surprisingly, this tran- 
sition is not shifted with respect to the pristine melt- 
ing transition. We interpret this transition as first-order 
melting of the vortices in the pristine bulk beneath the 
patterned surface, that results in a step in vortex density. 
The added vortex "tails" beneath the patterned surface 
force vortex "tips" through the surface, even in regions 
of reduced compressibility. We emphasize that this is an 
extremely unusual situation, in which the vortex tails lo- 
cated in the bulk are in a liquid phase, while their tips 
located near the surface are in a pinned, ordered phase, 
and therefore solid. 

The irreversibility line of the patterned regions is found 
to be shifted upward in H and T to above the pris- 
tine melting line T m , and displays step-like behavior, 
with almost no temperature dependence between match- 
ing fields B = nBff). These steps are consistent with a 
multi-vortex pinning scenario, with an estimated maxi- 
mum value of n max = 6 flux quanta pinned to each sur- 
face hole at H/Bcj, > 5 and T < 72 K. 

Applied shaking shifts the irreversibility line to lower 
field and temperature, enabling the observation of the 
first-order transition. However, no first-order transitions 
related to the surface holes, corresponding to depinning 
lines of interstitial vortices, or derealization of vortices 
pinned to holes, were observed. Further transport mea- 
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surements are necessary to determine if an additional de- 
localization line exists at higher temperatures, in anal- 
ogy to the derealization line of the vortices residing on 
columnar defects, that separates the interstitial liquid 
from a homogeneous liquicffl 
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